The NMSSM Solution to the Fine- Tuning Problem, Precision 
Electroweak Constraints and the Largest LEP Higgs Event Excess 
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We present an extended study of how the Next to Minimal Supersymmetric Model easily avoids 
fine-tuning in electroweak symmetry breaking for a SM-like light Higgs with mass in the vicinity of 
100 GeV, as beautifully consistent with precision electroweak data, while escaping LEP constraints 
due to the dominance oi h ^ aa decays with nia < 2mb so that a r~^r~ or jets. The residual 
^ 10% branching ratio for h ^ bb explains perfectly the well-known LEP excess at rrih 100 GeV. 
Details of model parameter correlations and requirements are discussed as a function tan f3. Compar- 
isons of fine-tuning in the NMSSM to that in the MSSM are presented. We also discuss fine-tuning 
associated with scenarios in which the a is essentially pure singlet, has mass nia > 30 GeV, and 
decays primarily to 77 leading to an h ^ aa ^ 47 Higgs signal. 



I. INTRODUCTION 

In the Standard Model (SM), electroweak symmetry 
breaking, whereby the W and Z bosons and the quarks 
and leptons acquire mass, gives rise to a Higgs boson, 
hsM- However, the value of rrihsM is quadratically sensi- 
tive to the cutoff scale of the theory. A, especially through 
top quark loops which give a one-loop correction of 
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where A is the high energy cutoff and = 176 GeV. For 
A of order the GUT scale, Mjj-, or the Planck scale, Mp, 
an extreme cancellation between the one-loop contribu- 
tion(s) and the bare Higgs mass is required in order that 
the physical Higgs mass be below a TeV, as required 
in order for the scattering of longitudinally polarized W 
bosons to obey unitarity in a perturbative fashion. 

Supersymmetric (SUSY) models, such as the Minimal 
Supersymmetric Model (MSSM), cure this naturalness / 
hierarchy problem associated with the quadratically di- 
vergent 1-loop corrections via the introduction of super- 
partners for each SM particle. Because the spin of the 
superpartners differs by 1/2 unit from that of the cor- 
responding SM particle, the 1-loop correction from the 
superpartner will cancel that of the SM particle once the 
energy scale being integrated over in the loop is above 
the mass of the (presumed to be heavier) superpartner. 
So long as the superpartners have mass somewhat below 
1 TeV (say ~ 500 GeV), the cancellation is not partic- 
ularly extreme and the hierarchy / naturalness problem 
associated with the quadratic divergences is ameliorated. 
However, there remains the question of how finely the 
GUT-scale parameters must be adjusted in order to get 
appropriate electroweak symmetry breaking, that is to 
say correctly predict the observed value of mz- It is here 
that LEP limits on a SM-like Higgs boson play a crucial 
role. 

Supersymmetric models most naturally predict that 
the lightest Higgs boson, generically /i, is SM-like and 



that it has a mass closely correlated to m^, typically ly- 
ing in the range < 105 GeV for stop masses < 500 GeV, 
with an upper bound, for example, of < 135 GeV in the 
MSSM for stop masses ~ 1 TeV and large stop mixing. 
If the stop masses are large, the predicted value of mz is 
very sensitive to the GUT scale parameters. Such sensi- 
tivity is termed 'fine tuning'. Models with minimal fine 
tuning provide a much more natural explanation of the Z 
mass than those with a high level of fine tuning. The de- 
gree of fine tuning required is thus quite closely related to 
the constraints on a SM-like /i, and these in turn depend 
on how it decays. 

The SM and the MSSM predict that h ^ bb decays 
are dominant and LEP has placed strong constraints on 
Zh Zbb. The limits on 
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are shown in Fig. 1 (from Ref. [1]). From this plot, one 
concludes that rrih < 114 GeV is excluded for a SM- 
like h that decays primarily to bb. In fact, because of 
the manner in which the analysis is done, at a first level 
of approximation this limit applies for an h that decays 
to any combination of 2b and 46. For B{h bb) ~ 
0.15 and B{h bbbb) ~ 0.8 (with r channels making 
up the rest) rrih ^ HO GeV is excluded. This will be 
important later. In the case of the CP-conserving MSSM, 
one always obtains B{h bb) > 0.88. For msusY ^ 
1 TeV, most of parameter space will yield rrih < 114 GeV 
and thus be ruled out by the SM-like Higgs LEP limit . 
The remaining part of MSSM parameter space either has 
at least one very large parameter, most typically a soft- 
SUSY-breaking stop mass close to a TeV at scale m^, 
or else large mixing in the stop sector. In the former 
case, one always finds that to predict the observed mz 
requires very careful adjustment, i.e. fine-tuning, of the 
GUT-scale parameters with accuracies better than 1%. 
In the latter case, fine-tuning can be reduced to the 3% 
level. To achieve small fine-tuning, let us say no worse 
than 10%, the soft-SUSY-breaking parameters that affect 
the Higgs sector should be well below a TeV, in which 
case the lightest CP-even MSSM Higgs boson would have 
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mass - 100 GeV. 

As suggested in [2] , the simplest way to allow a Higgs 
mass of order 100 GeV, thus making possible a light 
SUSY spectrum and low fine-tuning, is to modify Higgs 
decays so that the bb branching ratio is small and primary 
decays are to channel (s) to which LEP is less sensitive. 
This is very natural in models in which the Higgs sector 
is extended and Higgs to Higgs decays are kinematically 
allowed. The decay widths for Higgs to Higgs decays can 
easily exceed the very small width for the bb channel. 
The simplest supersymmetric model that gives rise to 
this possibility is the Next-to-Minimal Supersymmetric 
Model (NMSSM). The NMSSM yields a preferred value 
of rrih ~ 100 GeV purely on the basis of minimizing fine- 
tuning. A Higgs mass near 100 GeV is also strongly pre- 
ferred by precision electroweak measurements. Further, 
there is a well-known 2.3<j excess in the e+e~ Z -\- b's 
channel in the LEP data for M5/5 ~ 100 GeV when a 
final state that contains two or more b^s is assumed to 
contain exactly 2 b^s. If the Higgs decays only to bb then 
this excess and limits on the Z + b^s final state would 
apply to defined by 

Ce'^ = [glzJglzhsJB{h ^ bb) . (3) 

The excess is apparent in the higher observed vs. ex- 
pected limits for a test Higgs mass of rrih ~ 100 GeV 
shown in Fig. 1. This excess is particularly apparent in 
the l — CLb result (Fig. 7 of [1]) obtained after combining 
all four LEP experiments. 
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FIG. 1: Expected and observed 95% CL limits on C^^ from 
Ref. [1] are shown vs. ruh. Also plotted are the predictions 
for the NMSSM parameter cases discussed in [3] having fixed 
tan/3 = 10, Mi,2,3{mz) = 100,200,300 GeV that give fine- 
tuning measure F < 25 and < 2mb and that are con- 
sistent with Higgs constraints obtained using the preliminary 
LHWG analysis code [4]. 

In a previous paper [3] , we have shown that the above 
excess is consistent with a scenario in which the Higgs bo- 
son has SM-like ZZh coupling, but has reduced B{h 
bb) by virtue of the presence of h decays to a pair of 
lighter Higgs bosons, h aa^ where B{a ^ bb) is small. 



as is automatic if < 2m5 so that a t^t~ or light 
quarks and gluons. ^ (The importance of such decays 
was first emphasized in [5], and later in [6], followed 
by extensive work in [7-10].) For example, if the ZZh 
coupling is full SM strength, then ruh ~ 100 GeV with 
B{h bb) - 0.08 and B{h aa) - 0.9 fits the observed 
Z2b excess nicely. Meanwhile, there are no current lim- 
its on the Zh Zaa Zt^t~t^t~ final state for 
rrih ^ 87 GeV [11]. And limits in the case of a ^ jets 
run out at slightly lower mh- As already stressed and as 
described below in more detail, we are particularly led to 
the above interpretation of LEP data since fine-tuning 
within the NMSSM is absent for model parameters that 
yield precisely this kind of scenario [2, 3]. While vari- 
ous alternative interpretations of this excess in terms of 
a non-SM Higgs sector have been suggested [11, 12], the 
NMSSM scenario has the lowest fine tuning of any such 
scenario and has particularly strong theoretical motiva- 
tion. 

The NMSSM is an extremely attractive model [13]. 
First, it provides a very elegant solution to the /i prob- 
lem of the MSSM via the introduction of a singlet super- 
field S. For the simplest possible scale invariant form of 
the superpotential, the scalar component of S naturally 
acquires a vacuum expectation value of the order of the 
SUSY breaking scale, giving rise to a value of /i of order 
the electroweak scale. The NMSSM is the simplest su- 
persymmetric extension of the standard model in which 
the electroweak scale originates from the SUSY break- 
ing scale only. Hence, the NMSSM deserves very serious 
consideration. 

Apart from the usual quark and lepton Yukawa cou- 
plings, the scale invariant superpotential of the NMSSM 
is = XSHuHd + I i^S^ depending on two dimen- 
sionless couplings A, n beyond the MSSM. [Hatted (un- 
hatted) capital letters denote superfields (scalar super- 
field components).] The associated trilinear soft terms 
are XAxSHuHd + ^ k,A^S^ . The final two input pa- 
rameters are tan/5 = hu/hd and /ieff = As, where 
hu = {Hu)^ hd = {Hd) and s = {S). The Higgs sector 
of the NMSSM is thus described by the six parameters 
A , K , , A^, tan/3 , /ieff • In addition, values must be 
input for the gaugino masses and for the soft terms re- 
lated to the (third generation) squarks and sleptons that 
contribute to the radiative corrections in the Higgs sector 
and to the Higgs decay widths. 

The particle content of the NMSSM differs from the 
MSSM by the addition of one CP-even and one CP-odd 
state in the neutral Higgs sector (assuming CP conserva- 
tion), and one additional neutralino. The result is three 
CP-even Higgs bosons (/ii,2,3) two CP-odd Higgs bosons 
(ai,2) and a total of five neutralinos Xi,2,3,4,5- It will be 
convenient to denote the CP-even and CP-odd neutral 



■"^ If a ^ 66 is dominant, as occurs for rua > '2mi), then, as noted 
earlier in the text, rrih ^ HO GeV is required by LEP data [4]. 



3 



Higgs bosons of the MSSM as /i, H and A, respectively, 
while those of the NMSSM will be denoted by /ii, ^2 
and ai,a2, respectively. In the latter case, our focus 
will be on the lightest states hi and ai. The NMHDE- 
CAY program [14], which includes most LEP constraints, 
allows easy exploration of Higgs phenomenology in the 
NMSSM. 

In [2, 3], we presented a first study of the fine-tuning 
issues for the NMSSM vs. the MSSM. We define the 
fine-tuning measure to be 



MSSM: tanjS = 10, Mi,a.3=100.200.300 GeV 



F MaXpFp Maxp 



dlogmz 
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where the parameters p comprise all GUT-scale soft- 
SUSY-breaking parameters. 



II. COMPARISON OF THE MSSM WITH THE 
NMSSM 

In this section, we will consider scenarios associated 
with minimal fine tuning in the MSSM and the NMSSM. 
In the following section, we will give a broader overview 
of all types of NMSSM scenarios and will show how it is 
that one is lead to the NMSSM scenarios considered in 
this section. 

We discuss fine tuning for the MSSM first. In this 
case, the GUT scale parameters comprise: : Mi^2,3, 
B^, mfj, m|), m|, m|, mjj^, mjj^, At, A5, and A^. 
In principle, soft masses squared for the first two genera- 
tions should be included above, but they have negligible 
effect upon mz- In our approach, we choose m^-scale 
values for all the squark and slepton soft masses squared 
at scale m^, for the gaugino masses, Mi^2,3(^z), and for 
^t(^z), Ai){mz) and Ar{mz) (with no requirement of 
universality at the GUT scale). We also choose m^-scale 
values for tan/?, /i and m^; these uniquely determine 
B^{mz)- The vevs hu and hd at scale mz are fixed by 
tan/3andmz viam| = g'^{h1^h'^) (where = g'^^g''^). 
Finally, ^|/^(^z) and m?^^{mz) are determined from 
the two potential minimization conditions. [From here 
on, all parameters displayed without an explicit argu- 
ment are m^-scale values, although we sometimes give 
them an explicit (mz) argument for emphasis. All GUT- 
scale parameters will be specifically indicated using an 
explicit argument (Mf/).] We then evolve all parameters 
to the MSSM GUT scale (including ji and B^). Next, we 
shift each of the GUT-scale parameters in turn, evolve 
back down to scale mz-, and re- minimize the Higgs po- 
tential using the shifted values of /i, 5^, m|^^ and m\j^. 
This gives new values for hu and hd from which we com- 
pute a new value for mz (and tan/?). 

It is not difficult to understand why fine tuning is typ- 
ically large in the MSSM given LEP constraints. Mini- 
mization of the Higgs potential gives (at scale mz) 
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FIG. 2: Fine tuning vs. (top), rrih (middle) and At (bot- 
tom) for randomly generated MSSM parameter choices with 
tan;^ = 10 and Mi,2,3(mz) = 100, 200,300 GeV. Blue pluses 
correspond to parameter choices such that rrih < 114 GeV. 
Red crosses are points with rrih > 114 GeV. 



The m^-scale parameters can be deter- 

mined from the GUT-scale values of all SUSY-breaking 
parameters via the renormalization group equations. The 
result for tan^^ = 10 (similar to the tan/5 = 2.5 results 
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in Refs. [15, 16]) is 



rrir 



-2.{)i?{Mu) + 5.9M|(Mf/) + 0.8m|(Mt/) 
+0.6m^(Mf/) - l.2m]j^{Mu) 
-0.7M3(Mf/)A,(Mf/)+0.2A2(Mf/) + ... (6) 

All of the above terms aside from — 2/i^(M[/) and 
—2m?jj (Mu) arise from the RGE evolution result for 
2m|^ {Mu) — 2m\j Similarly, one can expand m^- 

scale values for soft-SUSY-breaking parameters in terms 
of GUT-scale parameters. In particular, one finds (at 
tan/3 10) 



At{mz) 
Ms{mz) 
m~{mz) 

In the above. 



-2.3Ms{Mu) ^0.24At{Mu) 
3Ms{Mu) 

5.0M|(Mt/) + 0.6mf (Mf/) 
-^0.2At{Mu)Ms{Mu) . 



(m~ + m~ ) 
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Unless there are large cancellations (fine-tuning), one 
would expect that 



(11) 



where mg± is similar in size to /i. We would need 
a very light gluino, and a rather light stop, to avoid 
fine-tuning. More precisely, if At{Mu) = 0, then it 
is clear from Eq. (6) that the minimum of F is deter- 
mined by the 5.9M^{Mu) term, which would give F ^ 6 
for Ms{Mij) = 300 GeV. Allowing for small positive 
At{Mu) reduces this minimum F somewhat, as will be 
illustrated below. Of course, in specific models you can 
also have correlations among the GUT-scale parameters 
that would reduce F. 

The problem is that the small At value required for 
minimal F does not yield a MSSM Higgs mass rrih above 
the 114 GeV LEP limit unless is very large (which 
causes a high level of fine-tuning, F > 175). To max- 
imize rrih at moderate m^, one should consider param- 
eters corresponding to \At/rnj\ ^ >/6, termed an 'm^- 
max scenario'. For such choices it is also possible to ob- 
tain ruh > 114 GeV. To simultaneously minimize F, the 
sign of At must be chosen negative. To understand this, 
we first note that rrih > 114 GeV can be achieved with 
-At > 500 GeV or At > 500 GeV and mj - 300 GeV. 
Given Eq. (9), this translates to -At{Mu) > 1 TeV or 
At{Mu) > 3 TeV, respectively. In both cases, F will be 
determined by the 0.2A^{Mu) term in Eq. (9), yielding 
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-At{Mu) > 1 TeV case, F > 30 is obtained. Similar val- 
ues of F can be obtained for much smaller \At{Mu)\ val- 
ues provided one allows for moderate fn~{Mu) < [17]. 

The above generic features are apparent in the numeri- 
cal results presented in Fig. 2 for the case of tan (3{mz) = 
10 and Mi,2,3(mz) = 100,200,300 GeV. We scan ran- 
domly over Atimz)-, Ai^imz)-, Ar{mz) and 3rd gen- 
eration squark and slepton soft masses-squared above 
(200 GeV)2, as weh as over \fi{mz)\ > 100 GeV, 
sign(/i) = zb and over ttia > 120 GeV. For such val- 
ues of rriA^ the h is quite SM-like and only allowed by 
LEP data if rrih ^114 GeV. If lower values of rriA are 
allowed, in particular rriA ~ 100 GeV, lower values of 
F ~ 16 can be obtained for experimentally allowed sce- 
narios. In these latter scenarios, the H is typically fairly 
SM-like but will have mass above 114 GeV while the h 
can have mass below 114 GeV by virtue having weak ZZ 
coupling. These scenarios are characterized by mixing 
among the Higgs bosons. Analogous mixed-Higgs sce- 
narios are also possible in the NMSSM. The MSSM and 
NMSSM mixed-Higgs scenarios will be considered in a 
separate paper [18]. The fine tuning in NMSSM mixed- 
Higgs scenarios have also been discussed in [19]. The 
main drawback of mixed-Higgs scenarios is that they re- 
quire adjustments in other parameters besides those nec- 
essary for correct electroweak symmetry breaking. 

Returning to Fig. 2, the top plot gives F as a function 
of m^. The latter enters into the computation of the 
radiative correction to the SM-like light Higgs mass rrih- 
In the middle plot, we display F as a function of rrih- 
And, in the bottom plot we display F as a function of 
At. We first of all note that the very smallest values 
of F are achieved for G [300 GeV, 400 GeV], rrih ^ 
90 - 105 GeV and At G [-400 GeV, 0]. As stated above, 
for rriA > 120 GeV, as considered here, the h is fairly SM- 
like in all its couplings to SM particles. Thus, points with 
rrih < 114 GeV, plotted as (blue) +'s, are excluded by 
LEP data, whereas those with rrih > 114 GeV, plotted 
as (red) x's, are not excluded by LEP. Although very 
modest values of F (of order F ~ 5) are possible for rrih < 
114 GeV, the smallest F value found for m/^ > 114 GeV 
is of order F ~ 34, as explained earlier. The increase of 
the smallest achievable F with rrih is illustrated in the 
middle plot. The modest F > 34 values are achieved for 
special parameter choices, namely, ~ 300 GeV and 
At ^ —500 GeV corresponding to a large ratio of At/rn^; 
see earlier discussion. A value of F ~ 34 corresponds to 
roughly 3% fine tuning. Generally speaking, however, it 
would obviously be nicer if the rrih ^ 100 GeV points 
with F ~ 5 were not excluded by LEP. 

We now compare these results to what is found in the 
NMSSM. Plots analogous to those for the MSSM appear 
in Fig. 3. Let us first first define our conventions (we 
follow those of Ref. [14]) and discuss a few theoretical 
points. The superpotential for the Higgs fields is 



Obviously the case of At{Mu) ^ 3 TeV case will corre- 
spond to very large fine tuning, roughly F > 180. For the 
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FIG. 3: For the NMSSM, we plot the fine-tuning measure 
F vs. m^, nihi and At for NMHDECAY-accepted scenarios 
with tan/3 = 10 and Mi,2,3(mz) = 100, 200, 300 GeV. Points 
marked by blue 'H-'s are consistent with LEP limits on the Z+ 
2b channel and the Z -^4b channel [11], considered separately, 
but not necessarily with LEP limits on the combined Z -\- 2b 
and Z -\-4b channels). Points marked by red 'x's escape LEP 
limits due to m^^ > 114 GeV. 
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For the soft SUSY breaking terms we take 
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(Above, we have not written the usual terms involving 
Higgs fields and quark/squark fields.) Assuming that the 
parameters of the potential are real, W and Vsoft together 
yield a full potential for the neutral components of the 
Hu, Hd and S scalar fields of the form 
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There are now three minimization conditions 
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which are to be solved for m'jj^ 



TTi^Hd ^1 terms 



of the vevs and other parameters appearing in Eq. (16). 
One combination of the minimization equations yields 
the MSSM-like expression for m| in terms of /i^, tan^, 
rn^Hu ^'hcl ^^^^ replaced by /ieff- However, a second 
combination gives an expression for /ieff in terms of mz 
and other Higgs potential parameters: 
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Eliminating /igff, we arrive at an equation of the form 
m%^Bm%^C = 0, with solution m| = -B±VB'^ - 4C, 
where B and C are given in terms of the soft SUSY break- 
ing parameters. A, and tan p. Only one of the solutions 
to the quadratic equation applies for any given set of pa- 
rameter choices. 

To explore fine tuning numerically, we proceed analo- 
gously to the manner described for the MSSM. At scale 
mz, we fixed tan^ and scanned over all allowed values 
of A (A > by convention and A < 0.7 is required for 
perturbativity up to the GUT scale) and a^, and over 
100 GeV < |/ieff| < 1.5 TeV, sign(/ieff) = ±. We also 
choose m^-scale values for the soft-SUSY-breaking pa- 
rameters Ax, A^, At = Ab, Ml, M2, M3, rriQ, m 



u-> 



m 



and m|;, all of which enter into the evolution equa- 
tions. We process each such choice through NMHDE- 
CAY to check that the scenario satisfies all theoretical 
and experimental constraints, with the exception that we 
plot some points that are consistent with LEP limits on 
the Z + 26 and Z + 46 channels considered separately as 
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in [11], but inconsistent with the LEP constraints on the 
Z -f h' s final states, where h' s = 2b-\- 4h. We shall return 
to this point shortly. For accepted cases, we then evolve 
to determine the GUT-scale values of all the above pa- 
rameters. The fine-tuning derivative for each parameter 
is determined by shifting the GUT-scale value for that 
parameter by a small amount, evolving all parameters 
back down to mz-, redetermining the potential minimum 
(which gives new values h'^ and h'^) and finally computing 
a new value for m| using m'^ = g'^iK^ + h'^'^^ 
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FIG. 4: For the NMSSM, we plot C^^ and as a functions 
of irihi for NMHDECAY-accepted scenarios with tan (3 — 10 
and Mi,2,3(mz) = 100,200,300 GeV. Point notation as in 
Fig. 3. ' ' 

Our basic results are displayed in Fig. 3. The density 
of points of a given type should not be taken as having 
any significance — it simply refiects the nature of the 
scanning procedures employed for the various cases. In 
particular, the scans used to obtain results presented in 
this section were designed to focus on parameter regions 
with F < 250. Further, we focused a lot of our scans on 
keeping only points with rria^ < 20 GeV. 

In Fig. 3, one sees a lot of similarity between the 
NMSSM plots and those for the MSSM, with differences 
to be noted below. Again, one finds that F < 10 {i.e. 
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FIG. 5: For the NMSSM, we plot the fine-tuning measure 
F vs. BR{hi aiai) and vs. rua^ for NMHDECAY- 
accepted scenarios with tan/3 = 10 and Mi,2,3(^z) = 
100, 200, 300 GeV. Point notation as in Fig. 3. 



no worse than 10% fine-tuning of GUT scale parame- 
ters) is easily achieved in the NMSSM for the present 
modest gluino mass of 300 GeV if the mean stop mass 
is of order 300 - 400 GeV which yields rrih, 100 GeV 
(for the case of tan^ = 10 — variation with tan^ will 
be noted later). The associated At values are of very 
modest size, lying in the range [—400 GeV, —100 GeV]. 
Further, as described in more detail later, the low-F sce- 
narios are once again such that the hi is quite SM-like 
as regards its couplings to lyVF, ZZ and //. The dif- 
ference between these plots and the earlier MSSM plots 
is that many (but not all, as we shall explain) of the 
rrih^ < 114 GeV points plotted escape all LEP limits. 
Of course, points with rrih^ > 114 GeV escape LEP lim- 
its simply by being above the maximum LEP-excluded 
mass. 

Let us next discuss in more detail the points with 
m/^^ < 114 GeV. They are plotted in Fig. 3 provided 
that they are consistent with LEP limits on the Z -\- 2b 
and Z ^ Ah channels, considered separately as plotted 
and tabulated in [11], but not necessarily the combined 
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Z + h's limit. All the plotted m/^, < 114 GeV points 
pass the Z + 26 and Z + 46 separate channel limits by 
virtue of the fact that B(h\ ^ a\a\) large enough that 
B(h\ bb) is sufficiently suppressed that Z-\-hi Z-\-bb 
lies below the LEP Z -\- 2b limit, while simultaneously 
Z -\- hi Z -\- aitti bbbb lies below the Z -\- Ab limit. 
The values of C^^ and C^^ are plotted in Fig. 4 as func- 
tions of rrihi. In the C^^ plot one finds two classes of 
points with < 114 GeV. The first class has large 
(but still below LEP limits for this individual chan- 
nel) by virtue of the fact that B{hi aiai) is large and 
rua^ > 2m}). The second class of points has zero C^^ 
since m^i < 2m5. F as a function of B{hi aiai) is 
shown in the top plot of Fig. 5. 

However, the plotted rrih^ < 114 GeV points with non- 
zero are mostly not consistent with the LEP data. 
NMHDECAY allows these points because it does not take 
into account the need to combine Z -\-2b and Z ^Ab final 
states in confronting the LEP limits, which are effectively 
(at least roughly) on the sum of these two final states. For 
those rrih^ < 114 GeV points with non-zero C^^, the sum 
of and C^^ is typically large and one expects such 
points to have too large a net Z + b' s rate, where b' s = 
2b + 46. Indeed, for the limited number of points that 
were analyzed using the full LEP Higgs working group 
code one indeed finds [4] that those points with rrih^ < 
110 GeV that have large B{hi aiai) and rria^ > 2m5 
(implying ai bb predominantly) are inconsistent with 
LEP limits on the net Z + b's rate. Without analyzing 
every one of our > 2m5 points using the full code, 
we cannot be sure that this same statement applies to all 
of them. 

With this proviso, we thus find that when m/^^ < 
114 GeV, one needs large B{hi aiai) and rua^ < 2m5 
to evade LEP limits. The bottom plot of Fig. 5 shows 
that it is easy to obtain very low F points that satisfy 
both criteria. (The less frequent occurrence of m^i > 
20 GeV points in this plot is purely an artifact of our 
scan procedure.) We will turn to a discussion of this in 
more detail shortly. 

As regards rrihi > 114 GeV points (which are not sub- 
ject to LEP limits), returning to Fig. 3, we find that F 
values as low as ~ 20 {i.e. only 5% tuning of GUT-scale 
parameters) can be achieved for the special choices of 
- 300 GeV and At - -500 GeV. This is the same 
region of stop parameter space that yields a minimum 
F - 34 with rrih, > 114 GeV in the MSSM. 

We have also found another type of point with low F 
and m/i^ ~ 100 GeV that escapes published LEP lim- 
its as follows. First, for these points B{hi aiai) 
is large, > 0.75, so that B{hi bb) < 0.2, imply- 
ing a perfectly acceptable LEP rate in the Zh Z2b 
channel. Second, the ai is highly singlet and decays 
mainly into two photons, B{ai 77) > 0.9. Thus, 
there is negligible contribution to the Zh Z4b chan- 
nel. Thirdly, m^^ is typically fairly substantial for these 
points, rria^ ^ 30 — 45 GeV. However, these points are 
highly fine-tuned in the sense that the highly singlet na- 



ture of the ai required for large B{ai 77) is very 
sensitive to GUT scale parameters. This is why they do 
not appear in the random scans discussed above. Locat- 
ing such points requires an extremely fine scan over a 
carefully chosen part of parameter space. We will give 
more details regarding these points later. 

We now briefly describe a third class of points that 
manage to have relatively low fine-tuning. Generically, in 
the NMSSM it is easy to have rrihi < 114 GeV without 
violating LEP limits simply by choosing parameters so 
that the hi has substantial singlet S component. In this 
way, the ZZhi coupling is suppressed and the e+e~ 
Z* Zhi production rate is reduced to an allowed level 
even if hi ^ bb decays are dominant. In such scenarios, 
it is typically the /12 that is the most SM-like CP-even 
Higgs boson, but > 114 GeV and LEP constraints 
do not apply to the /i2- We have performed a broad scan 
over NMSSM parameter space to look for and investigate 
the fine tuning associated with scenarios of this type. 
We find that not all the points of this type found in our 
scans are highly fine-tuned. There is a specific parameter 
region that produces points of this type that are only 
moderately fine-tuned and for which the hi escapes LEP 
limits by virtue of small ZZhi coupling. The lowest F 
value that we have found for such points is F ~ 16. In a 
separate paper [18], we will describe these scenarios and 
their fine-tuning in detail and compare to similar MSSM 
scenarios that are found when ttia ~ 100 GeV points are 
included in the MSSM parameter scans. 



III. THE LOW-FINE-TUNING REGION 

Armed with this overview, we now return to the pa- 
rameter region of the NMSSM that allows for the low- 
est possible fine tuning, as studied earlier in Section II 
for tan^ = 10. Here, we consider also tan/3 = 3 and 
tan^ = 50. These scans are focused very much on pa- 
rameter choices that can yield the lowest F values. The 
relevant plots are presented in Fig. 6. We present our re- 
sults in a somewhat different manner than in Sec. II so as 
to stress the remarkable preference for m/^^ ~ 100 GeV in 
order to achieve the very lowest F values at tan/3 = 10, 
with corresponding preferences for rrih^ ~ 101 GeV at 
tan^^ = 50 and rrihi ~ 90 GeV at tan^^ = 3. First, we 
stress that the above m/^^ values are the largest ones con- 
sistent with low F in an unbiased {i.e. before applying 
experimental constraints of any kind) scan over the part 
of parameter space that is simply theoretically consistent 
(see below). Once one imposes lower bounds on the stop 
(and chargino) masses, F shows a distinct minimum at 
the above nihi values. The preference for these values 
of rrihi to achieve low F becomes progressively more ap- 
parent as one imposes in addition: (a) LEP constraints 
on Higgs bosons, including the important Z -\-2b channel 
and Z -\- Ab channels, considered separately (as plotted 
and tabulated in [11]), but not the combined Z -\-2b and 
Z -\- 4b channels; and (b) LEP constraints on the com- 



8 



tan/?=10, Mi,2.3= 100,200,300 GeV 




.... I .... I .... I . 

I 60 60 100 

m^^ (GeV) 



tan/?=50, Mi,2.3= 100,200,300 GeV 



F 




40 60 60 100 



(GeV) 

FIG. 6: Fine tuning vs. ruhi for tan/3 = 10, tan/3 = 3 and 
tan/3 = 50 for points with F < 50, taking Mi,2,3{mz) = 
100, 200, 300 GeV. Small black x points are those obtained 
after requiring a global and local minimum, no Landau pole 
before Mu and a neutralino LSP. The green O's are those that 
in addition satisfy experimental limits on stops and charginos, 
but not necessarily Higgs limits. The blue D's are the points 
that remain after imposing all LEP single channel Higgs lim- 
its, in particular limits [11] on the Z -\-2b and Z -\-4:b channels 
considered separately. The yellow fancy crosses are the blue 
square points that remain after requiring ruai < 2mb, so that 
LEP limits on Z -\- b's, where b's = 2b -\- 46, are not violated. 



bined Z -\- 2b and Z -\- Ab channels. 

Let us first focus on the tan ^ = 10 case. Four different 
types of points are displayed. The black crosses show F 
as a function of rrih^ after requiring only that the sce- 
nario be theoretically consistent, but before any exper- 
imental constraints whatsoever are imposed. The most 
important components of the theoretical consistency are: 
(1) that the vacuum corresponds to a proper electroweak 
symmetry breaking vacuum at a true minimum of the 
potential; and (2) that the couplings remain perturba- 
tive during evolution up to the GUT scale. The black 
crosses already single out ^ 100 GeV as the point 
above which F rises rapidly. Black points with low ruh^ 
typically have a rather low value for m^^ that is clearly in- 
consistent with LEP and Tevatron limits. The minimum 
F for these low-m/^^ black cross points is fairly indepen- 
dent of rrihi . The (green) circles correspond to the black 
crosses that survive after imposing experimental limits 
on m^^ and m~± and similar non-Higgs constraints. We 
immediately see a striking preference for rrih^ ^ 100 GeV 
in order to achieve minimum F. The (blue) squares in- 
dicate the points that survive after requiring in addition 
that the scenario be consistent with LEP Higgs limits, 
including the Z -\-2b and Z -\- Ab final state limits consid- 
ered separately [11], but before imposing a limit on the 
combined Z + b's (b' s = 2b ^ 4b) final state. In the case 
of tan/3 = 10, these blue squares are the union of the 
F < 50 red x's and blue +'s of the middle plot of Fig. 3. 
The blue-square points now show a very strong prefer- 
ence for rrihi ^ 100 GeV, even before, but especially 
after, focusing on minimal F. The final large (yellow) 
crosses are the ruh^ < 114 GeV points among the (blue) 
square points that have m^^ < 2m5 so that there is no 
contribution to the Z -\- 4b channel from the hi aiai 
decay that, in turn, has a sufficiently large branching ra- 
tio to allow these points to escape the Z -\- 2b channel 
LEP limit. Now m/i^ ~ 100 GeV is clearly singled out. 

In the plots for the tan^^ = 3 and tan/3 = 50 cases, 
we did not bother to generate points with low mj^. So 
the black cross points in these cases simply indicate the 
presence of a few scenarios with m^^ above experimental 
limits but with m~± below existing limits or some other 
non-Higgs experimental inconsistency. The green-circle, 
blue-square and yellow-cross are as described above. 

The large number of blue-square points with very low 
F indicate that a significant fraction of the very lowest 
F scenarios are such that hi decays primarily into a pair 
of the lightest CP-odd Higgs bosons of the model, hi 
aitti. The yellow crosses show that low-F points with 
large B{hi aiai) and m^i < 2m5 are often found. 
For such points, ai t^t~ (or qq + gg if m^^ < 2mr) 
thereby allowing consistency with LEP constraints on the 
Z + b' s channel and, in many cases, the LEP excess in 
the hi bb channel for Higgs mass of order 100 GeV. 
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Is small niai natural? 



Given that low F can be easily achieved without vio- 
lating LEP constraints if m^i < 2m5, an important is- 
sue is whether obtaining small m^i requires fine-tuning 
of GUT-scale parameters. In fact, a light ai is natu- 
ral in the NMSSM in the A^.Ax limit. This can 
be understood as a consequence of a global U{1)r sym- 
metry of the scalar potential (in the limit A^^Ax 0) 
which is spontaneously broken by the vevs, resulting in 
a Nambu-Goldstone boson in the spectrum [6, 10, 20].^ 
This symmetry is explicitly broken by the trilinear soft 
terms so that for small A^^ Ax the lightest CP odd Higgs 
boson is naturally much lighter than other Higgs bosons. 
In fact, as discussed in depth in [10], the values of A^, Ax 
needed to have rua^ < 2m5, B{hi aiai) > 0.7 and low 
F are quite natural in the context of the NMSSM start- 
ing from small or zero values of A^{Mu) and Ax{Mu)- 
In particular, large B{hi aiai) is essentially auto- 
matic for typical RGE generated values of Ax and A^ 
{\A^{mz)\ < 10 GeV and \Ax{mz)\ < 200 GeV) and it 
is only a question of whether the requirement < 2m5 
is naturally achieved. In [10], we found that essentially 
no tuning of Ax and A^ is required in many model con- 
texts. For example, in the case of tan/3 = 10 tuning of 
the GUT-scale parameters needed to achieve appropriate 
Ax and A^ is likely to be minimal for scenarios in which 
the ai is about 10% non-singlet at the state-mixing, am- 
plitude level, i.e. 1% at the probability level. 

More precisely, let us define 



ai = cosOaAmssm + sin OaAs , 



(19) 



where Amssm is the usual two-doublet CP-odd state and 
As is the CP-odd state coming from the S field. Then, 
the mass of the lightest CP-odd Higgs boson in the sim- 
plest approximation is given by: 



3s 



3XAx cos^ 6 
2 sin 2p 



nA^ sin Oa 



(20) 



We see that the Ax contribution to rua^ is suppressed 
relative to the contribution for small cos Oa and large 
tan/5 and an appropriate balance between the contribu- 
tions is naturally achieved. In [10], we defined a measure 
called G that encapsulates the amount of tuning at the 
GUT scale that is likely to be needed to achieve small 
rria^. G is defined using 



Fa, 



_ Ax dml^ 



Fa. 



mi; 



dA, 



(21) 



as 



G = Min{Max[|F^J,|F^J],|F^, +F^J} . (22) 



^ The alternative for getting a light ai is to have a slightly broken 
Peccei-Quinn symmetry. However, the models with low F are 
not close to the Peccei-Quinn symmetry limit. 



As shown in [10], small G implies it is quite natural to 
get small even for fairly general Mf/-scale boundary 
conditions. For example, if Eq. (20) is approximately cor- 
rect, so that m^^ is linear in Ax and A^^ and if Ax and A^ 
are primarily sensitive to a single GUT-scale parameter 
p, then, if \Fax + Fa^ \ is small, sensitivity of m^^ to p is 
guaranteed to cancel. Nonetheless, the measure G should 
not be overemphasized since specific boundary condition 
choices can give small m^i even when G is large. 

In Fig. 7, we plot G as a function of cos Oa for param- 
eters choices that yield F < 15, taking tan^^ = 10, 3, 
and 50. The results are displayed using different shad- 
ings (colors) for different ranges of m^i, as delineated 
in the figure caption. For tan/? = 10, we see that G is 
minimal for |cos^yi| ~ 0.08, not much above the rough 
lower bound of \ cos Oa\ ^ 0.06. This lower bound is a 
direct consequence of the combined requirements that: 
1) B{hi aitti) be so large that B{hi bb) is small 
enough to escape the LEP limit on the Z -\- 2b channel; 
and 2) that < 2m5 so that ai does not decay to bb 
and there is no hi aiai bbbb decay contribution to 
the net Z-\-b^s channel. We also note that small G is only 
achieved for m^i > 2771^-, with m^^ > 7 GeV preferred. 

For tan/3 = 3, there is again a preference for larger 
rriai in order to achieve small G. However, small G can 
be achieved for a much larger range of cosOa- Of course, 
one should also notice that all the \ow-F solutions have 
cosOa < 0, with lower bound of | cos6'a| ^ 0.06. 

The bottom of the three plots shows that if tan /3 = 50 
then it is much more difficult to find solutions with low 
F that also have low G. The lower bound on |cos^a| 
needed to achieve large B{hi aiai) shifts downwards 
slightly to about 0.05. 

One should note that the coupling of the ai to bb 
is proportional to tan (3 cos Oa times the usual SM-like 
75 coupling strength. The lower limits on \ cos Oa\ at 
tan^ = 10, 3, 50 are such that \ tdiii /3 cos a\ ~ 1, ~ 0.6, 
~ 2.5. This means that the bb coupling is not particu- 
larly suppressed, and can even be enhanced with respect 
to the SM-like 75 value. (Of course, the ti coupling of the 
ai, proportional to cot P cos Oa^ is very suppressed.) The 
fact that the bb coupling is always significant implies that 
there is always a significant branching ratio for T jai 
(where the T can be the 15*, 2S or 3S state) so long as 
there is adequate phase space for the decay. The predic- 
tions for B{Tis jai) and further discussion appear in 
[21]. 



B. Dependence of F on the gluino mass 

The minimum value of F that can be achieved is, of 
course, dependent upon M3 (and is essentially indepen- 
dent of Ml and M2). Indeed, the largest GUT-scale pa- 
rameter derivative is very frequently that with respect 
to Ms{Mu). To explore this sensitivity, we have also 
performed a (somewhat less dense) parameter scan for 
the case of Mi,2,3 = 100,200,600 GeV at tan/3 = 10. 
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FIG. 7: For the F < 15 scenarios that are fully consistent with 
all LEP constraints, we plot G vs. cos^a taking Mi, 2, 3 = 
100,200,300 GeV and tan (3 = 10 (top), 3 (middle) and 50 
(bottom). The point coding is: black = 8.8 GeV < niai < 
dark grey (red) = 2mr < ruai < 7.5 GeV; light grey 
(green) = 7.5 GeV < rUai < 8.8 GeV; and darkest grey (blue) 
= TJlai < 2mr- 
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FIG. 8: Fine tuning vs. the Higgs mass for ran- 
domly generated NMSSM parameter choices except for fixed 
Mi,2,3(mz) = 100, 200, 600 GeV and a fixed value of tan (3 = 
10. Notation as in Fig. 6. 

The results for F as a function of m^^ are presented 
in Fig. 8. We find a minimum value of F ~ 20 at 
m/i^ ~ 104 GeV, the latter being somewhat higher 
than the ruh^ ^ 100 GeV location in the corresponding 
tan/3 = 10, Mi,2,3 = 100,200,300 GeV case. A SM-like 
hi with ruhi ^ 104 GeV is only consistent with LEP lim- 
its if B{hi aitti) is large and rria^ < 2m5 (the large 
(yellow) crosses). The rrih^ ^ 104 GeV location of the 
minimal F is less consistent with the M25 ~ 100 GeV ex- 
cess in the LEP data. However, to have m/^^ ~ 100 GeV 
in the Mi,2,3 = 100,200,600 GeV case is possible for 
F ~ 22, which is barely different from the F ~ 20 mini- 
mum value. 

One other point of interest is that < 2m5 points 
are more easily achieved at larger rrih^ ^114 GeV when 
Msimz) = 600 GeV than for Msimz) = 300 GeV. This 
can be understood by considering the special case of A = 
0.2, K = ±0.2 and tan/3 = 10. In this case, we find 

Ax{mz) 0.03A,{Mu) + 0.93Ax{Mu) - 0.35At{Mu) 

-0.03Mi(Mt/) - 0.37M2(Mt/) + OMMs{Mu), (23) 

A,{mz) - 0.90A,(Mf/)-0.11AA(Mf/) + 0.02A,(Mf/) 
+0.003Mi(Mt/) + 0.025M2(Mt/) - 0.017M3(Mf/) .(24) 

Let us consider the hi > case (by convention, A > 0), 
for which it can be shown that < is required [10] 
to get ml^ > 0. From Eqs. (24), (7) and (8), one 
finds that A^{mz) ^ OMMs{mz) + O.lA^(mz), im- 
plying that increasingly negative At{mz) is required to 
achieve < as M^iruz) increases. From Eq. (20), a 
small value of m^^ will be easily achieved in the present 
case of nA,^ < if Aximz) < so that the nA,^(mz) 
and \Ax{mz) terms tend to cancel. Eqs. (23), (7) and 
(8) imply Ax{mz) ^ -0.3Ms{mz) - l.5At{mz) which 
gives Ax{mz) < for increasingly negative At{mz) as 
Ms{mz) increases. In short, the larger M^{mz) is the 
more negative At{mz) can be while requiring small m^^. 
The more negative Atiruz)^ the larger stop mixing is at 
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fixed and therefore the larger rrihi • 



D. Properties of the heavier Higgs bosons for 
low-F scenarios 



C. Low-F scenarios and the LEP excess 

We will now discuss in more detail other proper- 
ties of the low-F scenarios with rrihi ^ 100 GeV and 

< 2m5, focusing first on the case of tan/3 = 10 and 
Mi,2,3(mz) = 100,200,300 GeV. First, we recah our 
earlier results from [3]. There, we studied in detail the 
F < 25 points from our earliest tan /3 = 10 scans as plot- 
ted in Fig. 1. The plot shows the C^^ predictions for 
all parameter choices in our scan that had F < 25 and 
m^i < 2m5 and that are consistent with the experimen- 
tal and theoretical constraints built into NMHDECAY 
as well as with limits from the preliminary LHWG full 
analysis code [4], which in particular incorporates limits 
on the Z + b's combined channel. Eight F < 10 points 
are singled out. As we have emphasized, these latter 
points cluster near rrih^ ~ 98^105 GeV. In [3], we found 
the remarkable result that not only are these F < 10 

< 2m5 points consistent with LEP limits, but also 
most are such that and B{hi bb) are appropriate 
for explaining the C^^ excess. We wish to emphasize that 
in our scan there are many, many points that satisfy all 
constraints and have m^^ < 2m5. The remarkable result 
is that those with F < 10 have a substantial probability 
that they predict the Higgs boson properties that would 
imply a LEP Zh ^ Z + 6's excess of the sort seen. We 
stress again that the F < 10 points with substan- 
tially above 2mb all predict a net Z + 6's signal that is 
ruled out at better than 99% CL by LEP data. Indeed, 
all such F < 25 points have a net hi 6's branching 
ratio, B{hi bb) -h B{hi aiai bbbb) > 0.85, which 
is too large for LEP consistency. In our larger scans, as 
represented by the C^^ results of Fig. 4, we see a huge 
number of m^^ < 2m5 points with approximately the 
correct C^^ to explain the LEP 100 GeV excess. 

For tan/3 = 50, Mi,2,3(mz) = 100,200,300 GeV, the 
preference for rrih^ ~ 101 GeV to achieve low F will 
again imply that many of the lowest F scenarios will 
provide a natural explanation of the Z -\- 2b LEP excess. 
At tan/3 = 3, the very lowest F values, F ~ 7 — 8, con- 
sistent with LEP limits are achieved for rrih^ ~ 95 GeV, 
as shown in the middle plot of Fig. 6. Such an hi mass is 
too low to provide a natural explanation of the Z -\-2b ex- 
cess. However, this same plot shows that the very slightly 
higher value of F ~ 10 is possible for rrih^ ^ 100 GeV. 
Thus, the LEP Z -\-2b excess is fully consistent with low 
fine-tuning scenarios that pass all LEP Higgs limits for 
all tan/3 > 3. (We have not explored still lower val- 
ues.) Of course, it is equally true that at tan^^ = 10 and 
tan/3 = 50, only a very modest increase in F would be 
needed for rrihi to take on a value that is not perfectly 
correlated with the location at M25 ~ 100 GeV of the 
Z + 26 LEP excess. 
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FIG. 9: Fine tuning vs. iriai (top) and vs. mf^+ (bottom) for 
points with F < 50 taking Mi,2,3(mz) = 100,200,300 GeV 
and tan f3 — 10. Point notation as in Fig. 6. 

Another interesting question is whether there is any 
correlation between F and or between F and the 
masses of the heavier Higgs bosons, m/^+, ma2^ ^/i2 and 
The plots for tan/3 = 10 appear in Figs. 9 and 
10. There, we observe that F depends very weakly 
on (once m^i < 27715), but that lower values of 

, , rrih2 •, are definitely preferred to obtain 
small F and that to obtain fully allowed yellow fancy 
crosses with F < 10 requires m/,+ G [300 GeV, 800 GeV], 
e [250 GeV, 750 GeV], rrih^ e [200 GeV, 600 GeV] 
and G [300 GeV, 800 GeV]. The corresponding plots 
for tan^ = 50 appear in Figs. 11 and 12. Again, we 
observe that F depends very weakly on mai, and that 
lower values of m/i+ , , ^/i2 5 ^/is are definitely pre- 
ferred to obtain small F. For tan/3 = 50, to obtain fully 
allowed yellow fancy crosses with F < 10 requires mh+ G 
[250 GeV, 750 GeV], m^^ e [250 GeV, 750 GeV], rrih^ e 
[200 GeV, 650 GeV] and G [250 GeV, 750 GeV]. 

These results have implications for the LHC and ILC. 
At the LHC, the main processes for producing and de- 
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FIG. 10: Fine tuning vs. 



(top), (middle) and 



TJih^ (bottom) for points with F < 50 taking Mi,2,3(^z) 
100, 200, 300 GeV and tan /3 = 10. Point notation as in Fig. 6. 



tecting these heavier Higgs bosons are gg ^ bbH (where 
H = a2^h2^hs) and gg bth~^ -\- gg bth~ , with, for 
example, H t^t~ . One finds [22] that detection be- 
comes possible when the hbH and bth'^ couplings are en- 
hanced by large tan/3. The mass ranges for the heavier 
Higgs bosons preferred for obtaining low F are such that 
if tan /3 = 10 they will be on the margin of detectability 
at the LHC, whereas if tan/3 = 50 they will certainly be 
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FIG. 11: Fine tuning vs. rriai (top) and vs. my^+ (bottom) for 
points with F < 50 taking Mi,2,3(mz) = 100,200,300 GeV 
and tan /3 = 50. Point notation as in Fig. 6. 



detectable. (At tan ;5 = 3 the small-F mass ranges for the 
^2, ^3, 0.2 are similar, but tan/3 is definitely too small for 
the above LHC modes to be detectable.) For the lowest 
part of the mass ranges, a signal for gg ^ bbH might also 
emerge at the Tevatron if tan ^ = 50. It is also important 
to note that the low-F mass ranges of the a2, /i2, ^3, 
are such that their pair production would mostly be out- 
side the kinematical reach of a = 500 GeV ILC, but 
that a substantial portion of the mass ranges are such 
that pair production would be possible at a = 1 TeV 
ILC. 



E. Features and parameter correlations for low-F 
scenarios 

We next turn to a detailed discussion of various cor- 
relations among the NMSSM parameters that are asso- 
ciated with low-F scenarios having large B{hi aiai) 
and rria^ < 2m5, i.e. the points indicated by the large yel- 
low crosses in the previous figures. We call such points 
"fully ok" . We first present some figures to illustrate how 
the fully ok points compare to points that are either ex- 
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FIG. 12: Fine tuning vs. 



(top), (middle) and 



TJih^ (bottom) for points with F < 50 taking Mi,2,3(^z) 
100, 200, 300 GeV and tan /3 = 50. Point notation as in Fig. 6. 



perimentally excluded or else have sufficiently large rrih^ 
(roughly rrih^ > 110 — 114 GeV) as to avoid LEP con- 
straints on the Z + b's channel. These latter points are 
the (blue) squares, (green) circles and black x's of the 
earlier plots. 

First consider Af. Plots of F as a function of At{mz) 
and At{Mu) are given in Figs. 13 and 14 for tan/? = 10 
and tan/3 = 50, respectively. These show that rather 



tan/?=10, Mi,2.3= 100,200,300 GeV 
:Vi i: r\^'\ ' — ' — i — 





A,{in^) (GeV) 
tan/?=10, Mi.2.3= 100,200,300 GeV 




A,(Mu) (GeV) 

FIG. 13: Fine tuning vs. At{mz) (top) and vs. At{Mu) 
(bottom) for points with F < 50 taking Mi,2,3(^z) = 
100, 200, 300 GeV and tan p = 10. Point notation as in Fig. 6. 



well-defined (and rather tan ^-independent) values are 
needed to achieve the very lowest F values, especially 
after imposing Higgs boson experimental limits. At scale 
mz the preferred At{mz) is of order —100 GeV. The 
corresponding At{Mu) is of order +600 GeV. The low- 
est F values are of course those associated with rrihi ~ 
100 GeV. This is consistent with our earlier discussion. 
The tan/3 = 10 points with large negative At{mz) val- 
ues that escape LEP limits by virtue of m/^^ > 114 GeV 
are the dark (blue) squares that begin at F ~ 20 and 
At{mz) - -500 GeV. 

In Figs. 15 and 16, we plot F as a function of /ieff 
(which in the case of the NMSSM is only defined at 
scale mz where EWSB has occurred) for the cases of 
tan^ = 10 and tan/3 = 50, respectively. As one could 
easily anticipate from Eq. (5) (with /i replaced by /ieff), 
fine-tuning is smallest for the the smallest values of /ieff- 
This figure also shows that |/ieff| values below about 
100 GeV are eliminated by the LEP limit on the mass of 
the lightest chargino. 

Next, let us examine the soft squark masses — Mq, 
Mf/, and Mjj — of the third generation. Values for these 
at scale mz are plotted in Figs. 17 and 18 for tan/3 = 10 
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FIG. 14: Fine tuning vs. At{mz) (top) and vs. At{Mu) 
(bottom) for points with F < 50 taking Mi,2,3(^z) = 
100, 200, 300 GeV and tan P = 50. Point notation as in Fig. 6. 



tan/?=10, M],2.3=100,200,300 GeV 




Meff(mz) (GeV) 

FIG. 15: Fine tuning vs. /Xeff for points with F < 50 tak- 
ing Mi,2,3(mz) = 100,200,300 GeV and tan/3 = 10. Point 
notation as in Fig. 6. 



and tan/3 = 50, respectively. We see that to obey Hmits 
on stop masses, there is a fairly definite lower bound on 
MQ{mz) and MD{mz), although low values for Mu{mz) 
are possible. And, to achieve F < 15 and satisfy all 
experimental limits requires all these soft masses to lie in 




Meff(mz) (GeV) 

FIG. 16: Fine tuning vs. /Xeff for points with F < 50 tak- 
ing Mi,2,3(mz) = 100,200,300 GeV and tan/3 = 50. Point 
notation as in Fig. 6. 

a very well defined band. The corresponding GUT-scale 
values are given in Figs. 19 and 20. Points with F < 15 
satisfying all limits again have soft masses squared at 
the GUT scale that fall within narrow bands (and are 
sometimes negative and sometimes positive). 

Next, we examine, for the case of tan p = 10, Ax and 
A^ at scale mz- Fig. 21 gives some results. The upper 
plot shows that F < 15 can be achieved for a wide range 
of Ax{mz)^ with points that obey all limits requiring a 
minimum value of \ Ax{mz)\ ^ 40 GeV. The middle plot 
shows that fully ok points require A^{mz) in a rather 
narrow band with 0.8 < \A^{mz)\ < 15. The bottom 
plot shows the correlation between A^{mz) and Ax{mz) 
that is required to get small m^i < 2m5, as discussed 
earlier. Note that either both must be negative or both 
positive for any point that is fully consistent with experi- 
mental limits. The lower bounds on their absolute values 
for the fully ok points — the large yellow crosses — are 
those required to have large enough B{hi — > aiai) to es- 
cape the Z + b's LEP limits for rrih^ ^ 100 GeV. Similar 
results are obtained for tan^^ = 50. All these results can 
be understood analytically as discussed in [10]. 

Fig. 22 shows the dependence on m^^ on Ax{mz) 
and A^{mz) in the case of tan^ = 10. One observes 
that large can be achieved for these same ranges of 
^a(^z) and A^{mz) just as easily as small ma^. It is 
just that cases with large rria^ > 2m5 and small F, which 
requires rrihi ~ 100 GeV, are not consistent with LEP 
limits on the net Z + b's channel, as we have discussed. 
Similar results are found for tan^ = 50. 

Plots of the GUT-scale parameters, Ax{Mu) and 
Af^{Mu)^ appear in Fig. 23. These show that the lowest- 
F scenarios that are fully consistent with experiment are 
often achieved for small values of these parameters. In 
terms of model building, these soft-SUSY-breaking pa- 
rameters are thus close to values associated with 'no- 
scale' soft-SUSY-breaking. 

Probably the most interesting parameter correlation 
is that regarding the soft-SUSY-breaking Higgs masses 
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FIG. 17: Fine tuning vs. MQ{mz), Mu{mz) and 
MD{mz) for points with F < 50 taking Mi,2,3(^z) = 
100, 200, 300 GeV and tan f3 = 10. Point notation as in Fig. 6. 
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FIG. 18: Fine tuning vs. MQ{mz), Mu{mz) and 
MoiTUz) for points with F < 50 taking Mi,2,3(^z) = 
100, 200, 300 GeV and tan f3 = 50. Point notation as in Fig. 6. 



squared at the GUT scale. These are plotted in Figs. 24 
and 25 for the cases of tan P = 10 and tan P = 50, respec- 
tively. These plots show that the fully ok scenarios with 
smallest F have very modest soft masses squared at the 
GUT scale, especially in the case of m'jj . Thus, some- 
thing close to a 'no-scale' model for soft Higgs masses 
squared at the GUT scale is preferred for low F. 

This preference for a 'no-scale' type of boundary condi- 



tion for the Higgs soft masses squared at the GUT-scale 
is further emphasized by the tan^^ = 10 plot of Fig. 26, 
where we overlap the values of m|^^(Mf/), m?^^{Mu) 
and m?g{Mu) for F < 50 yellow fancy cross scenarios 
of Fig. 24. 

We have said in many places that for the fully ok sce- 
narios the hi has quite SM-like h\ZZ coupling. This is 
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FIG. 19: Fine tuning vs. Mq{Mu), Mu(Mu) and 
Md{Mu) for points with F < 50 taking Mi,2,3(mz) = 
100, 200, 300 GeV and tan f3 = 10. Point notation as in Fig. 6. 
Negative values indicate cases for which Mq jj jj < in which 



case the plot gives —^J—Mq jj j^. 
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FIG. 20: Fine tuning vs. Mq{Mu), Mu{Mu) and 
Md{Mu) for points with F < 50 taking Mi,2,3(mz) = 
100, 200, 300 GeV and tan (3 = 50. Point notation as in Fig. 6. 



Negative values indicate cases for which Mq u jj < in which 



case the plot gives —'^—M^ ^ ^. 



illustrated in Fig. 27 where we plot F as a function of 



ghsmZZ 



(25) 



for the case of tan/? = 10 (results for tan (3 = 50 are 
similar). We see that the fully ok yellow fancy crosses all 



have \Cv\ ~ 1. In fact, for F < 50 scenarios, \Cv\ ~ 1 
also for the (blue) square points that are not also yellow 
fancy crosses, i.e. those points obtained if one only re- 
quires that the scenario is consistent with experimental 
limits that include the Z -\-2b channel (that is, before re- 
quiring rriai < 2mi) as needed to avoid the limits on the 



tan/?=10, Mi,g,3= 100,200,300 GeV 
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FIG. 21: F vs. Ax(mz) and At^{mz) (upper two plots) and 
^k(^z) vs. AA(mz) (lower plot) for points with F < 50 
taking Mi,2,3(mz) = 100, 200, 300 GeV and tan f3 = 10. Point 
notation as in Fig. 6. 



combined Z -\-b's channel). Suppressed \Cv\ values only 
appear in these plots if the Higgs experimental limits are 
removed. As discussed later, there are some very special 
points for which this is not true that will be considered 
in a follow-up paper. 
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FIG. 22: iTiai vs. Ax(mz) and A^(mz) for points with F < 
50 taking Mi,2,3(mz) = 100,200,300 GeV and tan/3 = 10. 
Point notation as in Fig. 6. 



F. Parameter correlations for the very lowest F 
points 

In this subsection, we consider at a still more detailed 
level the fully ok yellow points having large B{hi aiai) 
and rriaj^ < that also have F < 10. We will present 
results only for the case of tan P = 10 and Mi, M2, M3 = 
100, 200, 300 GeV. In general, the other tan/3 values give 
similar correlations aside from the shift in the value of 
ruhi that gives the lowest F value. In the plots presented 
in this section we will use blue +'s in place of the yellow 
crosses, since the latter do not display well on their own. 
Hopefully, there are few enough points on the following 
plots that the reader can match points from one plot to 
another. 

First, we show F as a function of rrih^ and m^i in 
Fig. 28. In the upper plot, we see again the preference 
for ruhi near 100 GeV. The lower plot shows that the 
very smallest F values occur at rria^ values above 2mr, 
implying that the ai t^t~ channel is the dominant 
ai decay. We note that the < rria^ part of the 

rria^ < 2m5 fully ok zone was also found in the companion 
paper [10] to be preferred in order to avoid fine-tuning 
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FIG. 23: F vs. Aa(Mc/) and A^{Mu) (upper two plots) and 
A^{Mu) vs. Aa(Mc/) (lower plot) for points with F < 50 
taking Mi,2,3(mz) = 100, 200, 300 GeV and tan (3 = 10. Point 
notation as in Fig. 6. 
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FIG. 24: Fine tuning vs. mH^{Mu), mji^{Mu) and 
m%{Mu) for points with F < 50 taking Mi,2,3(^z) = 
100, 200, 300 GeV and tan p = 10. Point notation as in Fig. 6. 



associated with getting small m^^ . 

One important prediction of any given parameter set 
is that for ^^(Z + h's), where 



dzZhsM 



B(hi bb) 



-^B{hi aiai) [B{ai bb)]'^ .(26) 

For the fully ok points one has < 2mb and thus 
^2(Z + b's) = C^^. More generally, ^^(Z + b's) is the net 
rate for LEP production of Z + 26 and Z -\-4:b final states 
relative to the rate that one would obtain for a SM Higgs 
boson which decayed entirely to bb. Of particular interest 
is the correlation between ^^(Z + b's) and F. Thus, an 
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FIG. 25: Fine tuning vs. mH^{Mu), m]j^{Mu) and 
m%{Mu) for points with F < 50 taking Mi,2,3(^z) = 
100, 200, 300 GeV and tan p = 50. Point notation as in Fig. 6. 



important comparison is the model prediction for + 
h' s) relative to the excess found at LEP in the vicinity of 
hh mass - 100 GeV. The value of ^^(Z + h' s) = as a 
function of rrih^ is plotted in Fig. 29. We observe that the 
points with F < 10 lie in the range C^^ ^0.2 with many 
of the very lowest F points having G [0.07,0.13], the 
range most consistent with the LEP excess at hh mass 
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FIG. 26: Fine tuning vs. the GUT-scale soft Higgs masses 
squared for points with F < 50 taking Mi,2,3(^z) = 
100, 200, 300 GeV and tan/3 = 10. Point notation: dark (red) 
im?jj^{Mu)] light grey (cyan) w?jj^{Mu)] darker grey (green) 
rris{Mu)- Points plotted are the yellow fancy cross points 
from Fig. 24. 
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FIG. 27: Fine tuning vs. the relative coupling strength, Cv = 
Qhiww / QhsM^w tan/3 = 10. Point notation as in Fig. 6. 



- 98 GeV. 

Next, we wish to illustrate the relative MSSM vs. sin- 
glet composition of the ai and hi for the F < 10 points. 
This composition has obvious implications for their cou- 
plings to SM particles. The more pure MSSM the hi is, 
the more SM-like will be its couplings. The more sin- 
glet the ai is, the more weakly it will be coupled to SM 
particles. In particular, its couplings to SM down-type 
fermions and leptons are given by tan /3 cos Oa times the 
SM-like weight in which hlh is replaced by hi^^h^ for ex- 
ample. These compositions are shown in Fig. 30. The 
upper plot illustrates that there is a lower bound on 
I cos^^l that arises from the joint requirements of large 
B{hi aitti) and m^^ < 2m5. As noted earlier, this 
guarantees that the hh (and r+r~) coupling strengths 
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FIG. 28: The upper plot shows fine tuning vs. nihi for the 
large (yellow) cross points (for clarity, we use blue +'s in 
their place in this and succeeding plots) with F < 10 taking 
Mi,2,3(mz) = 100,200,300 GeV and tan/3 = 10. The lower 

plot shows F vs. ITlai ■ 



of the ai are sufficiently large that the decays of the 
ai are dominated by the heaviest fermionic states, i.e. 
ai T^T~ for 2mr < rria^ < 2m5 and ai cc for 
2mc < rua^ < 2mr, with ai gg also being impor- 
tant. For < 2mc, ai gg is the dominant decay. 
Note the preference for cosOa ^ —0.1 for the very lowest 
F points. The lower plot of Fig. 30 shows the singlet 
component, sin 6^5-, of the hi for the fully ok solutions. 
The hi can be as much as 20% singlet at the amplitude 
level, but this means it is still 96% non-singlet in the 
amplitude-squared sense. As a result, all plotted points 
have \Cv\ ~ 1- The very lowest F points are clearly 
associated with very small sin ^5. 

W^e now turn to the GUT-scale parameters associated 
with F < 10 large yellow fancy cross points (plotted as 
blue +'s for these figures) that pass all experimental con- 
straints and the correlations among them. 

First, we consider n{Mu) and \{Mu) in Fig. 31. V^e 
see that the very lowest F values have fairly small \{Mu) 
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FIG. 29: The plot shows i^{Z + b's) (which equals Cl^ since 
Ceff = when rriai < 2mh) vs. F for the large (yellow) cross 
points with F < 10 taking Mi,2,3(mz) = 100,200,300 GeV 
and tan (3 — 10. 
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FIG. 30: The plots show F vs. cos^a (upper) and F vs. 
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FIG. 31: The upper plot shows fine tuning vs. X{Mu) for 
large (yellow) cross points with F < 10 taking Mi,2,3(^z) = 
100,200,300 GeV and tan/3 = 10. The lower plot shows 
K,{Mu) as a function of X{Mu). 



FIG. 32: The upper plot shows fine tuning vs. Ax{Mu) 
for fully ok points with F < 10 taking Mi,2,3(^z) = 
100,200,300 GeV and tan/3 = 10. The lower plot shows 
A^{Mu) as a function of Ax{Mu). 



and significantly larger n{Mu) ~ —0.4. 

W^e consider A^{Mu) and Ax{Mu) in Fig. 31. W^e see 
that the very lowest F values have fairly small Ax{Mu) 
and A,^{Mu)^ ie., as noted earlier, both are close to being 
consistent with no-scale boundary conditions at Mu- 

We consider mH^{Mu)^ rriHdiMu) and ms{Mu) in 
Fig. 33. We see that the very lowest F values have 
fairly small GUT-scale values for all the scalar Higgs mass 
squared values, again close to being consistent with no- 
scale boundary conditions at Mu- However, our scans 
did not locate any fully ok points for which mH^{Mu)-, 
rriHdiMu) and ms{Mu) were all simultaneously small. 
We are unsure at this time as to whether this is an arti- 
fact of limited computer time for scanning or something 
deeper. 



IV. MODERATELY LOW F POINTS WITH 
DOMINANT ai -f-f DECAYS 

Let us now turn to the special class of points mentioned 
in Section 11. These are the low-F points with a SM-like 
hi of mass ~ 100 GeV and for which > 30 GeV. For 
these points, B{hi aiai) > 0.75 and B{ai 77) > 
0.9. For this to occur, the ai must be highly singlet in na- 
ture so that the tree-level decays to fermion-antifermion 
are highly suppressed, in which case the chargino loop- 
induced decay to 77 can be dominant. (The relevant 
couplings are present even when the ai is purely singlet.) 
This combination of features allows consistency with the 
LEP limits on the Z-\-b's channel. We have not been able 
to determine if there are relevant limits on the Z + 47 
channel. This channel would have quite a high rate and 
most probably these relatively spectacular events would 
have been noticed. These points are also disfavored the- 
oretically since, as detailed shortly, a very high level of 
fine tuning of GUT-scale parameters is required in order 
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FIG. 33: The upper plot shows fine tuning vs. ruHui^u) 
for fully ok points with F < 10 taking Mi,2,3(^z) = 
100,200,300 GeV and tan/3 = 10. The middle plot shows 
tuh^ (Mu) as a function of tuhu {^u)- The bottom plot shows 
ms{Mu) as a function of mHu{Mu)- Our convention is that 
if an is negative then we plot —\^—M^. 



to achieve cos^^ ^10 ^ as required (for tan/3 = 10) for 
ai 27 to be the dominant ai decay channel. Nonethe- 



less, they should not be entirely discarded as a possible 
class and so we give some details regarding them in the 
following. These points were found using an extremely 
fine grid scanning approach of the type detailed in [10] 
with fixed msusY = —At = 300 GeV, /i = 150 GeV and 
tau;^ 10. 

Basic plots for this scenario appear in Fig. 34. The 
top plot of this figure gives F as a function of rrihi • We 
see again that rrih^ ~ 100 GeV gives the lowest F value, 
F ~ 11 in this case. The lower plot shows that m^i ^ 
30 TeV is required for this kind of scenario, with the 
lowest F obtained for ~ 30 GeV. Fig. 35 shows F 
vs. cosOa as the top plot, B{ai 77) as a function 
of cosOa as the middle plot and ^^(Z + b's) vs. rrihi as 
the bottom plot. The top plot is useful for correlating 
F with the value of cos^y^. However, note that there is 
some degeneracy: essentially the same values of F and 
cos Oa are sometimes obtained even though the basic scan 
parameters are different. The middle plot shows that 
B{ai 77) > 0.65 for these points, with B{ai 77) ~ 
0.65 for the F ~ 11 points. As expected, very small 
cosOa is required in order for the ai ^ 77 decays to be 
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FIG. 35: For the points with dominant ai ^ 77 decays and 
low we plot: F vs. cos^a (top); B{ai 77) vs. cos^a 
(middle); and ^^(ZH-6^s) vs. rrihi (bottom). Since F depends 
primarily on cos Oa, there are a number of points in the F vs. 
cos^A plot that are actually multiple repetitions of exactly 
the same F at a given cos^a value but with B{ai 77) and 
^^(Z + b's) varying slightly because of sensitivity to other 
parameters of the scan. 



dominant. The bottom plot of the figure shows ^^(Z + 
b's) as a function of m/^^. We observe that ^^(Z + b^s) is 
of the right general magnitude to explain the LEP excess 
for the lowest F points that have m^^ ~ 101 GeV. We 
note that ^^(Z + 6's) receives contributions from both the 
Z -\-2b final state from direct hi ^ bb decay and also the 
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FIG. 36: For the points with dominant ai ^ 77 decays and 
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a function of cos^a (bottom). (Many of the points have the 
same F^^^^ value.) 



Z -\- 4:b from hi aiai Ab where B{ai bb) < 0.35 
due to the competition from the ai 77 decays. 

The careful reader may wonder why it is that we can 
have small cosOa for these points whereas the m^^ < 
2m5 points have a lower bound on cos6>^. In fact, it 
is precisely the combination of the requirements that 
rriai < 2mi) and B{hi aiai) > 0.75 which forces a 
lower bound on cos^a- Values of cos^^i small enough 
to yield large B{ai 77) while at the same time 
B{hi aitti) > 0.7 is maintained are only possible for 
relatively large . 

The fine-tuning required in Ax and to achieve very 
small cos Oa can be quantified via the derivatives 



TT^COS 6 A 



d cos Oa Ax 
dAx cos Oa 



TTiCOSfc^A — 



() COS Oa A^ 
dA^ cos i 



(27) 

where all parameters are defined at scale mz- Under- 
standing of these quantities can be gleaned from the ap- 
proximate formula 



cos Oa 



Ml 



\v{Ax - 2ks) sin 2/3 
2\s{Ax + Ks) + inAf^s sin 2/3 
(28) 
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where we used ^ 

2Xs 



MgusY=-At=300 GeV, /x=150 GeV, tanj8=10 



sin 2(3 



{Ax + ^s) , 



= {Ax - 2ns) , 



2Xnv^ sin 2(3 + XA, 



v'^ sin 2/3 
^ 2s 



(29) 
(30) 

?>nA^s. (31) 



Eq. (28) shows that there wih be great sensitivity of 
cos6>^ to the value of Ax relative to 2i<is^ and al- 
most no sensitivity to A^. Both are confirmed by the 
numerical results we now present for fixed Mi^2,3 = 
100,200,300 GeV, tan/3 = 10, /ieff = 150 GeV,'At = 
-300 GeV, Ai, = Ar = 0, Mq^u,d,l,e = 300 GeV (for 
the relevant 3rd generation). Different points are ob- 
tained by scanning in X^ n^ Ax^ A^. (Obviously, many 
more 47 points could be found if the fixed parameters 
are allowed to vary. However, large < is essential to 
get small F for such points.) 

In Fig. 36, we plot Mf^/^li vs. cos (9a and 
vs. cosOa' The top plot shows that ^^2/^^! i^i^st 
be small for small cos^a (and there is a strong linear 
relation). The bottom plot shows that such small val- 
ues of cos^A imply rather large values of ^ . Given 
Eq. (23), high sensitivity to the m^-scale value of Ax 
implies a high level of fine-tuning for Ax (at scale ruz) 
with respect to Ax{Mu), At{Mu) and Ms{Mu). One 
should also note that the F^^^^ tuning measure for Ax 
is largest (~ —1500) for the point for which the EWSB 
fine tuning measure F is smallest and vice- versa. 

Finally, the two plots of Fig. 37 show that these points 
require largish tz and A that are fairly closely correlated, 
while A^ must be quite small. 



V. CONCLUSIONS 

There is strong motivation for a supersymmetric model 
with an extended Higgs sector containing one or more ex- 
tra Higgs singlet superfields. These motivations range 
from string theory model building, where it is known 
that SM-singlets are abundant in string theory compact- 
ifications, to the purely phenomenological, including the 
fact that adding anything other than SM singlets to the 
MSSM will typically destroy gauge coupling unification. 
In this paper, we have studied in detail the Next-to- 
Minimal Supersymmetric Model, which contains exactly 
one singlet Higgs superfield in addition to the two Higgs 
doublet superfields of the MSSM. We have shown that 
there is a portion of NMSSM parameter space with an 
abundance of attractive features, no outstanding prob- 
lems and which leads to an important set of predictions 
that should be taken quite seriously. There are many 



^ These mass squared matrix entries receive radiative corrections 
not shown here. 
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FIG. 37: For the points with dominant ai 77 decays and 
low F, we plot: k vs. A (top); and A^ vs. Ax (bottom). 



ways in which the NMSSM is a better benchmark the- 
ory than the MSSM, since it has important flexibilities 
that are currently leading to problematical issues for the 
MSSM. The attractive features of the NMSSM include: 

• a natural explanation for the /i parameter is pro- 
vided — since all superpotential couplings are di- 
mensionless in the NMSSM, the scale of /a is given 
by the scale of soft-SUSY-breaking, which (see be- 
low) can be well below a TeV; 

• the supersymmetric context provides a highly satis- 
factory solution of the naturalness / hierarchy prob- 
lem if the squark masses (in particular, the stop 
masses) and the gluino mass are well below a TeV 
(implying possible discovery at the Tevatron and 
very plentiful production at the LHC); 

• in particular, such squark masses imply that fine- 
tuning with respect to GUT-scale parameters is not 
required in order to obtain the observed value of 
mz and a light SM-like Higgs boson; 

• low squark masses imply that the lightest Higgs 
boson, the hi, will most naturally be SM-like in 
its couplings to SM particles and have a mass of 
order 100 GeV, close to the ideal value for satisfying 
precision electroweak constraints; 
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• LEP data is fully consistent with such an hi 
provided it decays mostly via hi aiai 
T^T~T^T~ (requiring 2mr < rUa^ < 27715) or 
hi aiai 4j (when rria^ < 2mr)^ where the 
ai is primarily the CP-odd component of the extra 
complex scalar Higgs singlet field; 

• an appropriately large value of B{hi aiai) is 
typically such that B{hi bb) ~ 0.1, thereby pro- 
viding a natural explanation for the event excess 
near 98 GeV in LEP data for the Z -\- bb channel; 

• it is quite natural for the ai to be lighter than 2m5 
while at the same time having sufficient hiaiai cou- 
pling for large B{hi aiai); 

• the optimal scenarios fit nicely with choices for the 
GUT-scale values of the soft-SUSY-breaking Higgs 
masses squared and A parameters that are quite 
modest in size, as might be associated with an ap- 
proximate no-scale model for SUSY breaking; 

• in the natural scenarios above, the heavier Higgs 
bosons of the model (two CP-even and one CP- 
odd neutral Higgs bosons and the charged Higgs 
boson) have relatively modest masses that would 
make them accessible at a hadron collider if tan p 
is large enough and mostly accessible at a 1 TeV 
linear e+e~ collider; 

• the light mostly-singlet ai must have a minimum 
coupling to the SM particles (through mixing with 
the non-singlet CP-odd state) that implies a lower 
bound, albeit small, on B{T 7^1); 

• the ai could allow for adequate annihilation in the 
early universe of very light neutralinos [23]. 

The attractiveness of this scenario suggests that the LEP 
groups should push a re-analysis of the Z4r channel in 
the hope of either ruling out the hi ^ aiai At sce- 
nario, or finding an excess consistent with it for rrih-^ in 
the vicinity of 100 GeV. Either a positive or negative 
result would have very important implications for Higgs 
searches at the Tevatron and LHC. We also stress that 
B factory experiments should attempt to search for a 
T jai signal down to the lowest possible branching 
ratio (the predicted minimum in the NMSSM context 
being of order 10~^). 

We speculate that similar results could emerge in other 
super symmetric models with a Higgs sector that, like the 
case of the NMSSM, is more complicated than that of the 
MSSM. Many such models can be constructed. Thus, 
much of the discussion above regarding Higgs discovery 
is quite generic. In general, there might be quite a few 



light a's, all of which could appear in the decay of a light 
SM-like h and all of which would provide potential signals 
in reanalyzed LEP data and in T ^ 7a decays. There 
is a potential gold-mine of discovery if one digs deeply 
enough. 

However, whether the a is truly the NMSSM CP-odd 
ai or just a lighter Higgs boson into which the SM-like 
h pair-decays, hadron collider detection of the h in its 
h ^ aa decay mode will be very challenging. Discovery 
modes that one can hope to demonstrate to be viable 
include: 

• WW fusion — WW ^h^aa^Ar; 

• tth production with h ^ aa ^ 4r; 

• diffractive production [24], pp pph^ with h 

aa 4r. This latter mode looks very promising 
[25]. 

Unfortunately, it seems very doubtful that viable discov- 
ery signals would be possible for the analogous modes 
with h ^ aa ^ A jet (that would be the only ones 
available if < 2mr). Although rria > 2mr is some- 
what preferred by naturalness arguments in the NMSSM 
case, one should be prepared for the possibility that the 
LHC will discover a plethora of supersymmetric parti- 
cles, and perhaps some heavy Higgs bosons (if tan/3 is 
large enough) but fail to see the SM-like light Higgs most 
closely associated with electroweak symmetry breaking. 
The only LHC evidence for its existence would then be 
that WW scattering would be found to be fully pertur- 
bative, as predicted if there is a light h with SM-like 
couplings to WW. 

At a linear collider, detection of e+e~ Zh produc- 
tion using the e+e~ ZX missing mass Mx approach 
will be completely straightforward. A 100 GeV h with 
SM coupling to ZZ will result in many events forming 
a sharp peak in Mx, quite independently of how the h 
decays. The decays can then be analyzed to see what is 
present and with what branching ratio. Detection of an 
h with unexpected decays at a photon collider will also 
be reasonably straightforward [26]. 
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